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ABSTRACT 

We follow the evolution of IC 10 X-l, the extragalactic binary hosting the most massive known stellar 
black hole with the mass of at least 23 M©. A massive companion of this black hole is a 35 M Q helium 
star that will very soon form another black hole. We demonstrate that this system will become a close 
double black hole binary with coalescence time of ~ 2 — 3 Gyr. We estimate that a detection rate of 
such systems is of the order of 0.5 yr -1 for initial LIGO and VIRGO. Additionally, we point out that 
the existence of a 23 M black hole in the small mctallicity environment of IC10 constrains the mass 
loss rates from massive stars to ~ 50% of currently used values. 

Subject headings: binaries: close — black hole physics — gravitational waves — stars: evolution — 
stars: neutron 

1. INTRODUCTION 

The interfcrometric gravitational wave observatories LIGO and VIRGO have already reached their design sensitivities 
and both are undergoing further improvements to reach better sensitivity. The most popular sources of gravitational 
waves that these experiments are looking for are coalescences of compact objects. Among these most attention has 
been paid to double neutron star systems (NS-NS). There is an observational evidence of their existence and their 
merger rates seem to warrant detection with the advanced interferometric experiments. The double black hole binaries 
(BH-BH) and black hole neutron star binaries (BH-NS) have received less attention. The direct detection of such 
systems in electromagnetic domain is difficult. From the theoretical point of view formation of such systems is not 
easy, as they have to pass through an unstable mass transfer phase which is not easy to survive for typical black hole 
masses of around 10M Q . 

The recent advances in X-ray instrumentation allow a study of X-ray binaries in the Local Group galaxies. IC10 X-l 
has already been discovered in ROSAT data by Brandt et al (1997). Bauer & Brandt (2004) have found an X-ray 
variability in a short Chandra observation. Clark & Crowthcr (2004) analyzed the possible optical counterparts of 
IC10 X-l(Crowther et al. 2003), and argued that it is a 35M Q WNE star. Subsequent longer Chandra observation 
lead to the discovery of X-ray periodicity (Prestwich et al. 2006). Prestwich et al. (2007) analyzed the X-ray and 
optical data of IC10 X-l, and found that it contains a black of a mass at least 24M in a binary with a w 35M Q 
companion. This result has been recently confirmed by Silverman & Filipenko (2008), who measured precisely the 
amplitude of the radial velocity of the companion. Thus black holes of stellar origin can reach much larger mass than 
previously thought, and moreover, they can be found in binaries with very massive companions. The other example 
is a binary M33 X-7 hosting a 16 M black hole with a 70 M Q companion (Orosz et al. 2007). We analyze the future 
binary evolution of IC10 X-l system using the Star Track binary evolution code. 

2. MODEL 

2.1. Population synthesis code 

Our population synthesis code, StarTrack, was initially developed to study double compact object mergers in the 
context of GRB progenitors (Belczynski, Bulik & Rudak 2002b) and gravitational-wave inspiral sources (Belczynski, 
Kalogera, & Bulik 2002a). In recent years StarTrack has undergone major updates and revisions in the physical 
treatment of various binary evolution phases, and especially the mass transfer phases. The new version has already 
been tested and calibrated against observations and detailed binary mass transfer calculations (Belczynski et al. 2008a) , 
and has been used in various applications (e.g., Belczynski & Taam 2004; Belczynski et al. 2004; Belczynski, Bulik & 
Ruiter 2005; Belczynski et al. 2006; Belczynski et al. 2007). The physics updates that are most important for compact 
object formation and evolution include: a full numerical approach for the orbital evolution due to tidal interactions, 
calibrated using high mass X-ray binaries and open cluster observations, a detailed treatment of mass transfer episodes 
fully calibrated against detailed calculations with a stellar evolution code, updated stellar winds for massive stars, and 
the latest determination of the natal kick velocity distribution for neutron stars (Hobbs et al. 2005). For helium star 
evolution, which is of a crucial importance for the formation of double neutron star binaries (e.g., Ivanova et al. 2003; 
Dewi & Pols 2003), we have applied a treatment matching closely the results of detailed evolutionary calculations. If 
the helium star fills its Roche lobe, the systems are examined for the potential development of a dynamical instability, 
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in which case they are evolved through a CE phase, otherwise a highly non-conservative mass transfer ensues. We 
treat CE events using the energy formalism (Webbink 1984), where the binding energy of the envelope is determined 
from the set of He star models calculated with the detailed evolutionary code by Ivanova et al. (2003). In case the CE 
is initiated by a star crossing the Hcrtzsprung gap (HG) we assume a merger and abort further binary evolution. This 
is due to the fact that there is no clear core-envelope boundary (and no entropy jump as for more evolved stars) in 
the interior structure of HG donors to facilitate the formation of a remnant binary system. As a consequence, a large 
decrease in the formation efficiency of close double compact binaries results (Belczynski et al. 2007). Most recently, 
the evolution of black hole spin was introduced to the model (Belczynski et al. 2008b) and the effect of massive twin 
binaries was investigated on double compact objects (Belczynski, Stanek & Fryer 2008c). For a detailed description 
of the revised code we refer the reader to Belczynski et al. (2008a). 

2.2. Host galaxy metallicity 

IC10 is a barred irregular galaxy in the Local Group at the distance between 600 and 800 kpc (Saha et al. 1996; 
Sakai, Madore & Freeman, 1999). It is undergoing a very rapid star formation and has a very high number of Wolf 
Rayet stars. IC10 has low metallicity: Lequeux et al. (1997) estimates it to be Z = 0.15Z©, but later studies by 
Massey et al. (2007) place it somewhere between the values for LMC and SMC. Thus we conservatively adopt the 
value ofZ = 0.3Z Q . 

3. RESULTS 

3.1. BH mass as a function of wind mass loss 

The mass of the BH in IC10 X-l is the largest of the known BHs of stellar origin. The measurement of such mass 
imposes interesting constraints on the stellar evolution. The mass of the stellar remnant depends on the metallicity 
and on the strength of the stellar winds. Here we use the formulas of Hurley, Pols, Tout (2000). They assume that 
the wind mass loss rate we use is 

M = f wmd (max(M R , M vw , M NJ , M WR ) + M LBV ) (1) 

where f W i„d is the wind scaling factor, and Mr is the mass loss rate in the giant branch (Kudritzki & Reimcrs 
1978), Mvw is the mass loss in the AGB phase (Vassiliadis & Wood 1993), Mnj is the massive star mass loss rate 
(Nieuwenhuizen & de Jager 1990), Mwr is the Wolf-Rayet star mass loss rate (Hamman & Kostoerke 1998), and 
Mlbv is the mass loss rate in the LBV phase (Humphreys & Davidson 1994). The mass loss from naked helium stars 
(Hurley, Pols & Tout 2000) is also scaled by the same factor f W i n d- We plot in Figure 1 the dependence of the mass 
of a black hole formed as a result of evolution of a single star with the initial mass 150 M Q as a function of the factor 
fwind for three values of the metallicity. The initial mass is close to the maximal mass of a star with a solar metallicity. 
The maximal mass of a BH strongly depends on the value of the wind scaling factor f W i n d and rather mildly on the 
metallicity. With the value of f w % n d = 0-5 the spectrum of black hole masses extends above the mass of IC10 X-l. 
This strong dependence has already been noted by Belczynski, Kalogera & Bulik (2002). 

3.2. The Future Evolution of IC10 X-l 

At present IC10 X-l has an orbital period 34.4 hr, black hole mass is estimated to be 23 — 33 M Q , while its companion 
is a helium star of a mass 17 — 35 M Q . The system is an eclipsing X-ray source with X-ray luminosity of 2 x 10 38 erg 
sec -1 an the lower limit on inclination was placed at 45 deg (Prestwich et al. 2007). 

For the edge on configuration (inclination: 90 deg) and the minimum donor mass of 17 M© one obtains minimum 
black hole mass: 23 M© (Prestwich et al. 2007). To form a black hole with such a high mass at the metallicity 
characteristic for IC10 (30% solar) it is required to reduce wind mass loss rates for massive stars. For standard wind 
mass loss rates a 150 M© forms a 13 M Q black hole, while for winds reduced by factor of 2 the same star forms a 23 M© 
black hole. Therefore, in our further considerations we will use the winds reduced by factor of 2 with respect to the 
fiducial value. To form a 17 M© helium star one requires a 50 M© star. The present black hole mass is probably very 
close to its formation mass, as the critical Eddington accretion rate limit for a 23 M Q black hole is 7 x 10~ 7 M Q yr _1 , 
while maximum potential accretion time is 1.5 Myr, allowing for accretion of 1 M©. The lifetime of a 150 Mq star is 
3.5 Myr, while a lifetime of a 50 M© star is 5 Myr. 

The evolution of the IC10 X-l from now on is very simple and is set by the present geometry of the system. The 
semi-major orbit axis is a = 18 R© (we have assumed a circular orbit), while radius of the massive helium star is 
^1 — 2 R© well within its Roche lobe Ri = 6.5 Rq. The lifetime of such a helium star is ~ 0.5 Myr and in this 
time the helium star loses ~ 5 M©. At the explosion helium star of ~ 12 M© and the CO core of mass ~ 9 M© goes 
through a direct collapse (Fryer 1999), i.e. no (or only very little) mass is lost and a massive 11 M© black hole forms. 
10% mass loss in neutrinos is assumed that induces a small eccentricity (e < 0.1) and slightly (< 10%) increase the 
orbit. The final double black hole binary with a 23 and an 11 M© black hole merges within ~ 2.2 Gyr, assuming that 
there was no natal kick involved in the formation of the second BH. 

However, should there be some asymmetry of the neutrino emission during the collapse or some mass ejection 
associated with black hole formation a moderate kick not exceeding 50km s -1 may be imparted on the newly formed 
BH. The orbital velocity in the pre-explosion binary is « 600km s -1 . Thus the system will not be disrupted by the 
explosion (Kalogera 1996), yet its orbit may become elliptical. This will influence the merger time of the system. We 
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have calculated the distribution of the merger times of the newly formed BH-BH binary assuming that the kick is 
drawn from the Maxwcllian distribution with a given width a v . For the assumed value of a v = 50km s -1 all systems 
remain bound, the median merger time becomes 2.3 Gyr, and for all but 1% it is smaller than lOGyrs. Increasing the 
width of the Maxwcllian to the a v — 100km s" 1 does not change the results significantly - the median merger time is 
2.5 Gyrs the probability of survival is 0.986, and 86% of the time the merger time is shorter than lOGyrs. Thus the 
possible kicks do not alter the future evolution of IC10 X-l significantly. 

3.3. Estimate of the coalescence rate 

The formation rate of IC10 X-l systems can be estimated as follows. We have seen so far one such system, and 
the Chandra or XMM sensitivity to detect them extends to distances within the Local Group. Let us assume that 
they could be detected up to the distance of r s , and only a fraction fl s of the sky has been searched for such systems 
therefore we can estimate the volume in which they are detectable as V s = f2 s 3 _1 r|? . The value of r s depends on the 
detectability of IC X-l like systems in X-rays, and also by the possibility of detailed spectroscopy of the companion 
in order to obtain the velocity curve. We estimate conservatively that it is about 2Mpc. We assume conservatively 
that almost all the sky has been searched for such binaries so il s = 4tt. The lifetime of the IC10 X-l in the X-ray 
bright phase is not longer than 1.5 Myr. This is a time in which accretion onto the first formed BH is allowed from its 
(non-BH) companion. The first BH is formed at 3.5 Myr, while the second BH is formed at 5 Myr after the formation 
of binary. This estimate was obtained under assumption that the companion mass is 17 M Q If companion is more 
massive, e.g., 35 M Q , the active X-ray phase is shorter than 1 Myr. We thus adopt the typical X-ray active lifetime 
of the order of t b s ~ 10 6 years. The formation rate of IC10 X-l like systems is then approximately estimated from 

p = v s -%l 

Currently LIGO can detect an NS-NS inspiral with the chirp mass of 1.2 M Q up to tens = 18Mpc. The chirp mass 
of the BH-BH system originating in IC10 X-l with the BH masses of 11 M Q and 23 M Q is M c hi rp ~ 13,6M Q and the 

5 /6 

detectability distance scales as the M c ^ i . Thus LIGO can detect mergers of BH-BH binaries originating in systems 

like IC10 X-l to the distance rsArs(M c ft,i rp /1.2) 5 / 6 = 172 Mpc. The formation rate of double black hole binaries in the 
volume observed by LIGO is therefore 
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For the Enhanced LIGO configuration the sampling distance for a NS-NS binary is 35Mpc, which brings the fiducial 
rate of equation [2] to 3.7yr _1 , while for the Advanced LIGO this distance will become 300Mpc, which implies the 
BH-BH merger rate of 2300yr~ 1 . 

Since the merger times for these binaries are below the Hubble time this is also the estimate of the merger rate 
detectable by LIGO and VIRGO. This estimate is rather conservative, since we assume that the star formation rate 
was constant over the past several billion years. Moreover the metallicity of galaxies was much lower a few billion 
years ago when the currently mcring systems were formed. This could also bring the current coalescence rate up. In 
the estimate of the rate we have neglected the fact that IC10 X-l is an eclipsing system. It is possible that there are 
additional such binaries in the Local Group yet they have not been analyzed because of lack of eclipses. Our estimate 
of the fiducial chirp mass of the system is also a lower limit obtained with the lower mass estimates both on a black 
hole and on its companion. Thus the above estimate must be considered as a rather conservative lower limit. On 
the other hand it is based on a single system and therefore one should treat it only as an indication of possibility of 
existence of a large population of merging massive black hole binaries. Moreover, in the estimate we have assumed that 
the properties of the Local Group is representative of a much larger volume. However the Local Group may contain 
an unusually large fraction of small metallicity galaxies which would bring the rate estimate down. The rate would 
also go down if the kicks received by the newly born massive black holes are much larger than currently believed, and 
reach the values typical for neutron stars. Thus lack of detection of such massive binary BH mergers by LIGO and 
VIRGO will imply that newly formed BHs receive large kicks, comparable to those of young neutron stars. This may 
be tested with the micro lensing experiments (Oslowski et al. 2008). 

4. DISCUSSION 

We have shown that the future evolution of the system IC10 X-l leads to formation of a double black hole binary 
with a merger time of the order of 2 — 3 Gyr. Such a system might be detectable by the intcrferomctric detectors like 
LIGO and VIRGO. Our estimate of the formation rate of such systems and the detection rate of their coalescences 
by LIGO is 0.5 yr -1 , which makes it very likely that such a merger can be found in the already gathered one year of 
the LIGO data. This rate exceeds even the estimated rate of NS-NS inspirals. The results do not depend on the kick 
that the newly formed BH may receive at its formation since the expected value of such kick is much smaller than the 
orbital velocity in the system. 

This estimate of the rate is, of course very crude. It relies on the properties of only one object. The rate calculation 
uses an assumption that the region sampled by Chandra is representative to the much larger volume in which such 
systems are detectable by LIGO or VIRGO. In particular we are assuming that the fraction of low metallicity galaxies 
is similar everywhere within the radius of ~ 170Mpc as in the Local Group. Yet this is the first estimate of the BH-BH 
merger rate based on observations of a binary which will become a double black hole within the next ~ 1 Myr. 
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The chirp mass of double black hole that forms out of IC10 X-l is very high ~ 20 M Q for a reasonable value of 
inclicnation. It is much higher than previously predicted for Galactic double black hole populations ~ 7 M Q (Bclczynski 
et al. 2007). This is due to the fact that the metallicity of IC10 is only 30% solar and we have decreased winds by 
factor of ~ 2 to allow for the formation of 24 M Q black hole observed in X-l. The double black hole binary that will 
form from IC10 X-l is very similar to the ones expected to form in globular clusters (Sadowski ct al. 2007). We stress 
that the inspiral templates used in search of gravitational radiation signals should be expanded to sample densely the 
region chirp masses close to and above 20 M©. 
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Fig. 1. — The mass of a black hole formed as a result of single evolution of a 150Mq star as a function of the wind scaling factor for 
three values of metallicity. 



